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Crystalline-Electric-Field Singlet-Triplet Configuration 
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Using the numerical renormalization group method, we investigate an extended Anderson 
model, in which correlated electrons with the ri (singlet)-^ (triplet) f 2 crystalline-electric- 
field (CEF) configuration hybridize with conduction electrons of ^(doublet) and T 8 (quar- 
tet) under cubic Oh symmetry, from a strong spin-orbit interaction limit. For the case of 
the parameters relevant to PrFe4Pi2, the system is under competition between the CEF 
singlet fixed point and the multichannel Kondo non-Fermi liquid fixed point arising from the 
quadrupolar coupling between the impurity with pseudospin 1 and the conduction electron 
with pseudospin 3/2 . We consider that this result reveals the origin of the heaviness of the 
effective mass and non- Fermi liquid behavior of the Pr-based filled skutterudite compounds. 

KEYWORDS: numerical renormalization group, f 2 singlet ground state, non-Fermi liquid, con- 
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Recently, a number of Pr-based filled skutterudite compounds have been investigated both 
experimentally and theoretically with much attention. Among them, PrOs4Sbi2 1 is attracting 
much attention because of time-reversal symmetry-breaking superconducting states 2 and mul- 
tiple phase diagrams. 3 Since the crystalline-electric-field (CEF) level for PrOs4Sbi2 has been 
well established 4-6 by both neutron scattering and ultrasonic experiments, it is important to 
know what happens under this CEF level scheme, in which the ground state is a T\ singlet 
with the lowest excited state is a T4 triplet, f-electron physics under CEF exhibits many 
interesting properties such as a two-channel (quadrupolar) Kondo effect 7 ' 8 and competition 
between the CEF and Kondo- Yosida singlet. 9 

In this Letter, we study an extended impurity Anderson model with the f 2 -CEF scheme of 
T\ — hybridizing with conduction electrons of Tj and T% symmetry. Recently, several studies 
using only a Tj conduction electron band have been carried out, 10,11 since it appears according 
to band calculations 12 that the conduction electrons in Pr-based filled skutterudites mainly 
come from the pnictogen molecular orbital with Ai u (i.e., Tj) symmetry. However, except for 
PrRu4Sbi2, there are certain amounts of contributions from the transition metal atoms, 13 e.g., 
PrOs4Sbi2 and PrFe4Pi2- It is natural to think that each filled skutterudite compound can be 
distinguished not only by the CEF levels but also by the variety in the conduction electron 
bands. Then it is not unrealistic to investigate a model including both and Tg conduction 
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electrons particularly as an impurity problem. 

First, let us start by introducing an extended impurity Anderson Hamiltonian appropriate 
to the present case. We restrict ourselves within the Hilbert space that includes T\ (singlet) 
and T4 (triplet) states for the f^-CEF configurations and T7 (doublet) and Tg (quartet) states 
for f^CEF ones under Oh symmetry. Hybridizations t>r 7 and t>r 8 are expressed easily if we 
work in a j-j coupling scheme with J = 5/2 manifold. Although there is no mixing between 
T4 and T§ states in Oh symmetry (the finite mixing is a striking feature in Th symmetry 14 ), 
we expect that such a obstacle will not hinder us from discussing the qualitative properties of 
a system with small deviations from Oh symmetry such as PrFe4Pi2- Even in the j-j coupling 
scheme, which is valid in a strong spin-orbit limit, we can include the effects of Hund's rule 
coupling in our model by considering the f 2 states with the total angular momentum J = 4. 

Thus, we state the model Hamiltonian as 

H = H c + H{1 + Hp + -ffhyb, (1) 



E[E e " 7c l c ^ + E e ' 8 «l«4 (2) 
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+ E <>Ul flr n(f 2 ri| + h.c.}, (5) 

M=±§,±± 

where the conduction electrons of Tg^y) are expressed as fermions dkuickv) with pseudospin 
3/2 (1/2) and the same representation for the f-electron operator a^(c u ) is used, rj is the CEF 
splitting between the T\ and T4 states. 

We can express two relevant f 2 states with the total angular momentum J = 4 (ri, T4) 
as the eigenstates of the total pseudospin S tot : Ti — > S tot = and T 4 — ► S 1 * * = 1 (see Table 
I). Then, we can reduce significantly the calculation time and maintain a high accuracy in 
the numerical renormalization group (NRG) calculation. 17 Namely, we can take advantage of 
the existence of three conserved quantities: the total pseudospin S tot , the total charge Q tot 
relative to the half filling, and an additional quantity, total helicity h tot , which is a quantity 
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Table I. A set of J = | (f 1 ) and J = 4 (f 2 ) states used for present NRG calculations, v =|, J. and 
/i = ±i, ±|. 



similar to the angular momentum in C3 V , C4 V symmetry, 15 and is denned as 

h tot =mod{^(aJ ( a jU + ^4 /1 a A;/1 ),2}. 



(6) 



The helicity is introduced to compensate for the lack of charge conservation law for each 
orbital. It is noted that h tot does not include c\c v . Then, f^-Fi states with J = 4, consisting 
of the states etc] and (ala_i — 4^3)1 f° r example, have a definite helicity h tot = 0. We 

2 _ 2 2 ~2 

can use h tot for sorting states in the NRG calculation, enabling us to diagonalize the block 
Hamiltonian, at each NRG step, faster by about four times. 

In this Letter, we use the discretization parameter A = 3 and perform calculations by 
retaining up to 1000 states, 17 while we have checked that the results obtained by retaining 
2000 states give no difference. For simplicity, we set the density of states of conduction electrons 
to be constant and the bandwidth 2D centered at the Fermi energy for both and Tg orbitals. 
For a later purpose, we define the hybridizations v? 7 and vr s as 



v Tr = -2v r 2 



Other hybridizations in eq. (5) are related to vp 7 and vp s according to Table I. Throughout 
this paper, energy is measured in units of D. 

Next, to make the physics of this extended Anderson model (EAM) clear, we transform 



eq. (1) to a Kondo-like Hamiltonian in the usual way, 18 assuming e f J 
2 = 1,4 and j = 7, 8. The result is 

-^Kondo = [</47S7 + J^8] ' S4 + JqQ8 ■ Q4 

+H pot + H n 

1 



Evi > K 7(8) | for 



1 



'47 



\VT 7 \ 



+ Hf2, 
2 



J48 



(8) 
(9) 



where 2 Jq = J^s- S7, ss and S4 represent the pseudospin of the conduction electron with T7, Fs 
symmetry and that of the impurity with T4 symmetry, respectively, qs and Q4 represent the 
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Fig. 1. T-dependence of 5i m p, entropy due to an impurity, for EAM (open symbols) and KM (filled 
triangles). Parameters used for EAM are -Er 7(s) = —0.9, Ep 1 — — f .5, vr 8 = 0.3 and vr 7 = 0.0. 
Although the result of EAM is obtained by setting vr 7 = 0, the qualitative features are the same as 
in the case of Vr 7 ^ 0. For KM, we use eq. (8) without H { 2 regarding J 48 and Jq ( J48 = Jq = 0.2) 
as independent parameters. 

pseudo-quadrupole of the T§ conduction electron and that of the impurity for T4 symmetry, 
respectively. It is noted that qs and Q4 are represented by the corresponding pseudospin 
operators as follows: {S x S y + S y S x , S y S z + S z S y , S Z S X + S X S Z , S% - S y , (3S Z - S 2 )/V3}. H pot 
in eq. (8) represents the potential scattering at the impurity site. We ignore the k-dependence 
of all parameters and the small renormalizations of the f 2 -level Er 4 and Er 1 ■ -ff m ix m eq. (8) 

contains terms such as ~ a\c^ (r^~ | representing the scattering that mixes I\ and and 

2 

will turn out to play an important role in the stabilization of the CEF singlet. The essential 
difference from the LS scheme 10 ' 11 is that it is not possible to rewrite the H m [ x term using the 
"X" operator 11 in the present model. This is because there is no triplet state constructed by 
two Tf electrons in the present scheme. It is noted that J47 < , eq. (9), which means that 
this coupling is irrelevant as in an usual Kondo model (KM). This point is consistent with 
the LS scheme. 

If we ignore H m ^ and Hp, we can obtain an impurity spin 1 and a conduction electron 
spin 3/2 multichannel KM with quadrupolar interactions, which was discussed by Koga et 
al. If we retain only the J48 term in eq. (8), we can construct an exact boundary conformal 
field theory (BCFT), 19-21 and the non- Fermi liquid (NFL) fixed point spectrum is in complete 
agreement with the NRG result. 16 If both J48 and Jq terms in eq. (8) were considered (no 
particle- hole symmetry), another NFL fixed point would be obtained, as shown by Koga et al 
in their NRG calculation. We call this NFL fixed point "S" following Koga et al. They argued 
that at the "S" fixed point, there exists a pseudo SU(3) symmetry according to the result 
of the two-loop renormalization group analysis. However, the NRG energy spectrum seems 
to have no SU(3) symmetry. Detailed analysis of the "S" fixed point will be investigated 
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elsewhere. 22 When we include all the terms in eq. (8), it is easier to treat the original EAM of 
eq. (1) by NRG. Taking the above knowledge into account, it is expected that the interplay 
between the CEF singlet and the NFL is important in the present model. 

We show the result of Si mp , the entropy due to an impurity, in Fig. 1 together with that 
of KM without H mix and Hp. It is noted that S imp for rj = 0.001, 0.005, and 0.01 in EAM 
(1) , and that in KM (8) merge into the same value ~ ks log 2 in the low temperature limit, 
suggesting that the ground state of EAM with those parameters is the T4 triplet state. As 
seen in Fig. 2, this is confirmed by the T-dependence of the number of doubly occupied states 
Ni and N±, which are defined as 



In Fig. 2, the T-dependence of Si mp is shown together with those of N\ and N4 for typical 
cases with the CEF singlet and the triplet ground state. A small deviation in Aq and N4 from 
those ideal values is due to the effect of hybridization with conduction electrons. The meaning 
of liniT^o <Si mp ~ ks log 2, which is not exactly equal to ks log 2, is unclear for the moment. 
We can also see the same NRG spectrum in both KM and EAM at the "S" fixed point. The 
NRG spectrum confirms that fee log 2 does not arise from the contribution of a simple spin 



The "phase diagram" in the T — r] plane is shown in Fig. 3, the phase boundary of which is 
obtained by estimating the crossover temperature of Si mp . Although we draw a phase diagram 
only for the parameter set vr 7 = in Fig. 3, the qualitative features for vr 7 7^ remain the 
same. The overall features are summarized as follows: a) The NFL "S" fixed point is possible 
in the case T s S> rj, where T s is the characteristic temperature corresponding to the energy 
gain in the T4 triplet state relative to the Ti level due to the Kondo effect, b) The CEF singlet 
fixed point is possible in the case T s <C rj. There exists the critical value rj c that separates a) 
and b). All the fixed points are explained by the shifts in the f 2 devels. 

Even in the region of the CEF fixed point, it is interesting to see whether the quasiparticles 
become heavy or not. We estimate Ci mp /T, specific heat divided by the temperature due to an 
impurity around the CEF singlet fixed point as shown in Fig. 4. It is a common feature that 
there exists a Schottky peak in C\ mp /T that corresponds to releasing S\ mp as /cb log 4 — > 0, 
and Ci mp /r approaches a saturated value as T — > 0. Ci mp /T at T — ► becomes large as t>r g 
increases (Fig. 4(a)) and as vp 7 decreases (Fig. 4(b)). These are precisely a reminiscence of 
the "S" fixed point (i.e., the effect of the critical point between the CEF singlet and the NFL 
ground state in Fig. 3). It is apparent from Fig. 4 that the heaviness of the quasiparticles is 
due to the hybridization with T$ conduction electrons. On the other hand, t>r 7 works to help 
in stabilizing the CEF singlet ground state, i.e., it renormalizes rj to become large. 



|f 2 r 1 )(f 2 r 1 



(10) 




(ii) 



1/2. 
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Fig. 2. Occupation number of f^-state together with the entropy due to an impurity. The upper 
panel is for the CEF singlet fixed point (vr 8 = 0.3) and the lower one is for the S NFL fixed point 
(vr 8 = 0.35). Other parameters are £r 7(8) = -0.2, E Tl = -1.5, r] = 0.0005 and v Tl = 0.3. The 
same symbols are used in both panels. The high-temperature entropy is not ~ fee log 4, because 
of the truncation procedure. The correct value 5i m p = kBlog4 is recovered when we keep 2000 
states in the NRG calculation. 
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Fig. 3. "Phase diagram" in T — r\ plane for v-p 8 — 0.37, v-p 7 = 0,-Epi = —1-4 and Ep 7(s) — 0.2. The 
critical value for r\ is determined as r\ c — 0.02 for these parameters. At r\ = rj c , the levels of the 
f 2 singlet and triplet states are interchanged. The crossover temperatures between quartet and 
triplet states, and between quartet and singlet states are proportional to |ry — r/ c \ as indicated in 
the figure. 



The parameter set of a small v r g or a large CEF splitting ij corresponds to the PrRu4Sbi2 
case, and the other limit corresponds to the PrFe4Pi2 case. The heaviness of the quasiparticle 
of the latter compound is related to the "S" fixed point. At this "S" fixed point, the scaling 
dimension of the leading irrelevant operator given by BCFT is 1/6 so that Ci mp /T shows 
a divergent behavior as C- mip /T ~ T~ 2 / 3 (see Fig. 5). 16 ' 22 However, in the experimentally 
attainable temperature region, we observe the — logT -like dependence of C mip /T, as seen in 
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Fig. 4. (a) Specific heat divided by temperature, Ci mp /T, vs T for Ti ground state, for a series 
of the hybridization vr 8 with other parameters fixed as «r 7 = 0.4, r\ = 0.0001, Er 1 = —1.5 
and Er 7(s) = —0.9. r\ is renormalized to ~ 0.005 for these parameters, (b) That for a series of 
the hybridization vr 7 , with other parameters fixed as vr 8 = 0.3, rj = 0.0001, Er 1 = —1.5 and 
Er 7(8) = -0.9. If we set D/k B = 10 4 K, C imp /T ~ 0.57 x (value of ordinate) mJ/mol • K 2 in both 
(a) and (b). 



the inset of Fig. 5. In this temperature region, our results are similar to the prediction based 
on the two-channel Kondo model, even though we did not assume that the ground state is 
the f 2 -F23 doublet. It is noted that the origin of the heaviness of the quasiparticles is different 
from the LS scheme. 11 

Finally, let us make a comment about the effect of symmetry. In this paper, we have 
restricted ourselves to Oh symmetry, because of certain technical reasons. The f 2 -Fs state with 
J = 4 is represented in our pseudospin language as 



•'21 & + 2 Z ' 2 V ZL 2 2 



T° : -±=[c\J 3 - c\a\] + */H[ t x0 t 3 _ a t a t i]; (13) 
V42 1 2 1 2 V 21 - - ~~ 



2 2 2 



^■.^[UaU + fclaK]-^aU 3 . (14) 

~'21 Z, 1 2 Z 1 2 V Zl 2 ~2 

It is noted that states (12)-(14) are not the eigenstates of S tot , 5 lot and h tot . If we include 
these states, considering Th symmetry (under which T4 and T§ mix), the "S" fixed point 
should become unstable because anisotropy exists there in the exchange interaction, breaking 
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Fig. 5. Specific heat divided by temperature, Ci mp /T, for T4 ground state. The parameters used are 
-^t 7(8) = —0.2, E Tl — —1.5, v rs — 0.5, v r? — 0.3 and rj — 0.0005. The result is consistent with the 
prediction of BCFT, C imp /T ~ T" 2 / 3 in the limit T -> 0. Inset: C imp /T around T ~ 10~ 4 . It is 
noted that this region is experimentally attainable and not extremely close to the "S" fixed point. 

the SU(2) symmetry. Such features of an unstable NFL fixed point can be seen in a recent 
NRG study reported in ref. 23, in which they discussed a case including T3 and T5 states in 
addition to the ground state quartet (Fx and F^) as the low-energy f 2 states (thus breaking 
the SU(2) symmetry in the present paper). However, the precise identification of the unstable 
NFL was not performed in the study detailed in ref. 23, although it seems to be related to the 
"S" fixed point discussed in the present paper. Another crucial point is that there is no triplet 
state constructed by two Tj electrons even if we work under symmetry. For this reason, 
the Kondo effect caused by Fj conduction electrons 11 is suppressed in the case of the strong 
spin-orbit interaction we are considering here. From these facts, we can expect that in the case 
of a strong Hund's rule coupling, it is possible for the Kondo effect of T-? conduction electrons 
to occur; on the other hand, in the case of strong spin-orbit interactions, the criticality related 
to the "S" fixed point creates strong renormalizations, if contributions from the Tg conduction 
electrons are important. 
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